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SWME Coiiaboration 

We report on recent progress in the calculation of Bk using HYP-smeared staggered fermions on 
the MILC asqtad lattices. Our main focus is on the continuum extrapolation, which is done using 
(up to) four different lattice spacings — a « 0.12, 0.09, 0.06 and 0.045 fm. Since Lattice 2010, 
we have reduced the statistical errors on the a « 0.09 fm lattices by a factor of ^ 3, and roughly 
doubled the size of the a w 0.045 fm ensemble. We find that these improvements have a very 
significant impact on the continuum extrapolation, with the a « 0.12 fm data lying outside the 
range of applicability of simple functional forms. Hence we use only the three smallest lattice 
spacings to perform the extrapolation, finding Ba: = Bk(RGI) = 0.725 ±0.004(stat) ±0.038 (sys). 
This value is consistent with our published value from 2010 (based the three coarsest lattice 
spacings), but has smaller errors. 
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1. Introduction 

The calculation of the kaon mixing matrix element Bk is one of the successes of lattice QCD. 
Results with all errors controlled are available with different fermion discretizations, the most 
recent entry being that using Wilson fermions [|l|]. For a summary, see Refs. and [^]. There 
is some tension between these results, which needs to be resolved in order to know whether the 
Standard Model can describe CP-violation in kaon mixing. 

In this proceedings we update our results for Bk- These are obtained using improved staggered 
fermions, specifically HYP-smeared valence quarks on asqtad sea-quarks. We describe here results 
obtained with chiral fitting functions from SU(2) staggered chiral perturbation theory (SChPT), 
which give our most reliable results [Q]. We compare our results to those obtained a year ago and 
presented at Lattice 2010 [^]. In companion proceedings we update the analysis based on SU(3) 
SChPT ^ and discuss strategies for dealing with correlations in chiral fits [^. 

Table |l] lists all the ensembles on which we have calculated Bk, and notes which results have 
changed in the last year. In particular, since Lattice 2010 we have accumulated higher statistics on 
the C2, C5, Fl, Ul ensembles and added a new measurement on S2. Of these, the most important 
updates are those on Fl and Ul, since they are used (together with C3 and SI) to do our continuum 
extrapolation. We focus on this issue here. The other updates, as well as the yet unanalyzed 
ensemble S2, provide information on the sea-quark mass dependence at multiple lattice spacings. 



a (fm) 


ami /cinis 


geometry 


ID 


ens X meas 


Bk{^i = 2 GeV) 


status 


0.12 


0.03/0.05 


20^ X 64 


CI 


564 X 1 


0.556(14) 


old 


0.12 


0.02/0.05 


20^ X 64 


C2 


486x9 


0.568(16) 


update 


0.12 


0.01/0.05 


20^ X 64 


C3 


671 X 9 


0.566(5) 


old 


0.12 


0.01/0.05 


28^ X 64 


C3-2 


275 X 8 


0.570(4) 


old 


0.12 


0.007/0.05 


20^ X 64 


C4 


651 X 10 


0.563(5) 


old 


0.12 


0.005/0.05 


24^ X 64 


C5 


509x9 


0.566(4) 


update 


0.09 


0.0062/0.031 


28^ X 96 


Fl 


995 x9 


0.527(4) 


update 


0.09 


0.0031/0.031 


40^ X 96 


F2 


850 X 1 


0.551(9) 


update 


0.06 


0.0036/0.018 


48^ X 144 


SI 


744x2 


0.537(7) 


old 


0.06 


0.0025/0.018 


56^ X 144 


S2 


198 X 9 


-NA- 


new 


0.045 


0.0028/0.014 


64^ X 192 


Ul 


705 X 1 


0.530(7) 


update 



Table 1: MILC asqtad ensembles used to calculate Bk- ami and am^ are the masses, in lattice units, of the 
light and strange sea quarks, respectively, "ens" indicates the number of configurations on which "meas" 
measurements are made. The results for Ba:(NDR,2 GeV) are obtained using the SU(2) SChPT fits (4X3Y- 
NNLO) discussed in the text. 

2. Chiral fits 

In our numerical study, our lattice kaons are composed of valence (anti)quarks with masses rrij^ 
and my- On each MILC ensemble, we use 10 valence masses: 

amx,amy = am^ X n/\0 with «= 1,2,3, 10, (2.1) 
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Figure 1: X-fits to BK{l/a) for the Fl ensemble, with 1 measurement/configuration (left) and with 9 mea- 
surements/configuration (right). Xp is the mass-squared of the xx valence pion, and we fix amy — 0.031 
(the heaviest value). In the right panel, the blue curve includes finite-volume (FV) corrections, while the 
red curve does not. FV corrections are not included in the left panel. Extrapolated results (diamonds) have 
known taste-breaking discretization errors removed and thus do not lie on the curves. 



where ams is the strange sea quark mass. In our standard fits we extrapolate to arn^ using the 
lowest 4 values for amx (the "X-fit" — done at fixed anty), and then extrapolate to mf^^^ using the 
highest 3 values of aruy ("Y-fit"). As described in Ref. [Q], these choices keep us in the regime 
where we expect next-to-leading order (NLO) SU(2) ChPT to be reasonably accurate. The X-fits 
described here are done to the form predicted by NLO partially quenched SChPT (and given in 
Ref. 10]), augmented by a single analytic term of next-to-next-to-leading order (NNLO). The amy 
dependence (which is not controlled by ChPT) is very close to linear and we use a linear fit for our 
central values. We dub this entire fitting procedure the "4X3Y-NNLO fit", and use it for our central 
values. Other choices (e.g. NLO vs. NNLO) are used to estimate fitting systematics. 

Examples of the X-fits are shown in Fig. |l]. Here we compare results on the Fl ensemble from 
Lattice 2010 with those from our present dataset. In each panel, a fit to the SChPT form is shown, 
together with the result after extrapolating am^ — > cim^^^^ and removing known discretization errors. 
For details of this procedure, see Ref. [Q]. 

The addition of 8 extra measurements per configuration reduces the statistical errors by 3, 
as expected. In addition, the central values shift down by about la, leading to a correspondingly 
lower value of Bk- We have also now incorporated finite volume corrections, as predicted by NLO 
SChPT, into the fitting function 



3. Continuum Extrapolation 

We do our continuum extrapolation using the four lattices with ami /arris = 0.2. Figure^ shows 
how this extrapolation is changed by our updated results. Clearly, the most significant change 
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Figure 2: Continuum extrapolation of Ba'(M — 2 GeV) versus (in units of 0.01 fm^). The left panel 
shows results from Lattice 2010, with linear and quadratic fits. The right panel shows our updated results, 
with constant and linear fits to the leftmost three points. 
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Figure 3: Continuum extrapolations of = 2 GeV) with different fitting choices. Fits are described in 
the text. 



is for the Fl ensemble, where the downward shift rules out the possibility of a straightforward 
extrapolation using all four points. We thus first consider extrapolations based only on the smallest 
three values of a. The figure shows results both for a constant and a linear fit, which have good 
X^/d.o.f as listed in Table 

We now describe our attempts to fit all four points to a reasonable functional form. As ex- 
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fit type 


#data 


fit function 


;^^/d.o.f 


;^^/d.o.f. (Bayes) 


const 


3 






0.70 


-NA- 


lin3 


3 


ci + C2a^ 




0.73 


-NA- 


lin4 


4 


ci + C2a^ 




6.35 


-NA- 


quad 


4 


ci + C2a^ + csd^ 




1.44 


4.39 


a2a2g2 


4 


ci +C2a^ + C3a^ai(a) 




1.66 


4.16 


a2g4 


4 






3.49 


4.48 


a2a2g2g4a4 


4 






-NA- 


3.86 



Table 2: Functional forms used for continuum extrapolations. Fits to 3 data points exclude the largest value 
of a. Bayesian fits and are described in the text. 

plained in Ref. the expected dependence is a linear combination of a^, a^as{l/a) and aj{l/a), 
together with higher order terms. Given that we have only four points, we began by using three- 
parameter fitting functions, in particular those labeled quad, a2g4 and a2a2g2 in Table |2[ in ad- 
dition to the two-parameter linear fit 1 in 4. Examples of the resulting fits are shown in Fig. |3[ The 
left panel shows unconstrained fits to all points using the linear (lin4-blue curve) and quadratic 
(quad-brown curve) functions, as well as reproducing the constant and linear fits to the three 
smallest values of a from Fig. ^ Both the lin4 and quad fits are problematic. The former has 
a poor (see Table ^, although the value of the coefficient C2 ~ (0.36 GeV)^ is reasonable. The 
quadratic fit has a reasonable but has unphysically large values for the coefficients. In par- 
ticular, we expect \c2\ = Aj and |c3| = A3 with A2 ~ A3 « Aqcd ~ 0.3 GeV, while the fit gives 
A2 ~ 0.49 GeV and A3 0.95 GeV. The latter value is unreasonably large. Similar problems occur 
for the a2a2g2 and a2g4 fit forms. 

Thus we have repeated the a2a2g2, a2g4 and quad fits imposing Bayesian constraints. We 
augment the x'^ in the usual way, with the expected central values of the coefficients C2 and C3 being 
zero, while the expected standard deviations are 

= 2A^ for a2a2g2, a2g4 & quad, (3.1) 
Oc, = (2A^), 2 & (2A**) for a2a2g2, a2g4 & quad, (3.2) 

with A = 300 MeV. The resulting fits are shown in Fig. ^ (b), with the augmented ;^^/d.o.f given in 
Table ^ The red curve shows the a2a2g2 fit, the green curve the a2g4 fit and the blue curve the 
quad fit. For reference we also show the unconstrained lin4 fit in brown. We find that all three 
constrained fits are poor, with the a2a2g2 and quad fits differing little from the lin4 fit. 

We have also attempted a fit to ci + 020^ + c^a^ as{a) + C4as{a)^ + €50"^ , which is possible as 
long as one uses Bayesian constraints on C2-5. The resulting fit, shown by the black curve in Fig. || 
(b), remains poor, with X'^/d.o.f. given in Table ^. Thus we conclude that we cannot describe the 
data from all four values of a with any of our fit functions if we insist on physically reasonable 
coefficients. Most likely this indicates that more terms are needed in the fit functions. 

Fortunately, the uncertainty in the correct global fit form is not that important for the value 
in the continuum limit. We take the result from the const fit for our central value, and use the 
difference with the Bayesian a2a2g2 fit for an extrapolation systematic. As Fig. || (b) shows, we 
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would obtain essentially the same systematic error were we to use either of the other two Bayesian 
fits or the lin4 fit. 

4. Updated Result and Error Budget 

Our updated result for Bk from the SU(2)-SChPT analysis is 

BA- = Bi?(RGI) = 0.725 ±0.004(stat)±0.038(sys). (4.1) 

Here we use the 4X3Y-NNLO fit to valence quark mass dependence and the const fit for the 
continuum extrapolation. 



cause 


error (%) 


memo 


status 


statistics 


0.58 


4X3Y-NNLO fit + const 


update 


matching factor 


4.4 


(Ul) 


il 


discretization 


2.7 


diff. of const and a2a2g2 (Bayesian) 


update 


fitting (1) 


0.92 


X-fit (C3) 


i 


fitting (2) 


0.08 


Y-fit (C3) 


i 


ami extrap 


0.06 


diff. of (C3) and linear extrap 


i 


GMs extrap 


0.5 


diff. of constant vs linear extrap 


i 


finite volume 


0.59 


diff. of y = oo and FV fits 


update If] 




0.14 


r\ error propagation (C3) 


i 


h 


0.38 


132 MeV vs. 124.4 MeV 


e 



Table 3: Error budget for Bk{2 GeV) obtained using SU(2) SChPT fitting. The "memo" column gives a 
brief description of the error estimate, while "status" indicates whether the error is updated or the same as in 
Lattice 2010 [§] or our paper [|]. 

The error budget for this result is given in Table ||. Most of the errors are estimated as in 
Ref. 10], and many are unchanged from that work, as indicated in the "status" column. The major 
changes in the last year are in the statistical and discretization errors. The former has been sub- 
stantially reduced compared to the 1.4% quoted at Lattice 2010 |^. By contrast, the discretization 
error has substantially increased (from 0.1%), because of our poorer understanding of the contin- 
uum extrapolation. Previously we used the difference between the result on the Ul ensemble and 
the continuum value as an estimate of this error, while here we use the difference between the 
const and a2a2g2 fits. 

We have also changed our method of calculating finite volume errors, but this has a minor 
impact on the final error. 

Combining the statistical and systematic errors in quadrature, our updated result has a 5.3% 
error This is slightly larger than the 4.8% error in our result from Lattice 2010 (Bk = 0.720 it 
0.010 ±0.033 [j5|]). The increase is due to the larger discretization systematic, which overwhelms 
the reduction in the statistical error. Still, the overall error is changed very little, because our 
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dominant systematic remains the truncation error introduced by our use of 1-Ioop perturbative 
operator matching.' 

In summary, improving the statistical errors has brought to light a difficulty in the continuum 
extrapolation that was previously masked. Although we have not fully understood the depen- 
dence, this has relatively little impact on our final result, because the extrapolation is anchored by 
the result from the smallest lattice spacing, a result which has not changed significantly. Never- 
theless, we intend to try more elaborate fits to improve our understanding of the dependence. 
We also need to revisit the sea-quark mass dependence, since the lower value for Bk on the Fl 
ensemble is now significantly different from that on the F2 ensemble. This is in contrast to the 
very weak sea-quark dependence on the coarse ensembles. Our most important need, however, is 
to reduce the truncation error, which we aim to do both using non-perturbative renormalization and 
two-loop matching. 
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